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In a series of papers, de Vega and Sanchez claimed that 
the thermodynamic hmit of a self-gravitating system can 
be taken by letting the number of particles, iV, and the 
volume, V. tend to infinity keeping the ratio N/V^^^ 
constant [jj]. This limit, which I call diluted following 
the terminology of the first paper of , is different from 
the usual thermodynamic limit, where the density N/V 
is kept constant. The relevant variable for the diluted 
limit, which can be found by naive dimensional analysis, 
is = Gm^N/V^/^T. 

Recently, I proved rigorously that the diluted limit 
does not give a well defined thermodynamic limit 0: 
the relevant thermodynamic potentials are not extensive 
and the thermodynamic quantities suffer from the same 
problems as in the usual thermodynamic limit. For in- 
stance, the free energy scales with N^/'^. However, in 
"The Cluster Expansion for the Self-Gravitating gas and 
the Thermodynamic Limit" de Vega and Sanchez 
continue arguing that the diluted limit gives extensive 
thermodynamic potentials and well behaved thermody- 
namic quantities at sufficiently high temperature, i.e. for 
•q > rjc, where rjc depends on the thermodynamic en- 
semble as well as on the geometry of the system bound- 
ary. To mantain these statements, these authors try to 
show that "the statements made in ref.[7] [reference [31 of 
the present paper] have crucial failures which invalidate 
the conclusions given in ref.[7]". However, the argument 
they give is based on a misunderstundingof the proof of 
nonexistence of the diluted limit given in '2| , and is easily 
refuted, as will be seen in the following. 

Let us briefly remenber the proof of nonexistence of 
the diluted limit given in P|. Consider a system of N 
classical particles, enclosed on a region of linear size R 
(so that V = B?) interacting via a gravitational poten- 
tial conveniently regularized at short distances. In the 
diluted limit we have N — > oo, i? — > cx), with N/R con- 
stant. The variable 77 of ref. Q is, by definition. 



inequalities, it is proven in that 



V = 



Gm^N 
RT 
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since V — R^ \s the volume of the region available for 
the system. Now, let us consider a region of linear size 
i?0) with N ^ R^, enclosed in the available space of the 
system. Note that Rq <^ R. Using a simple sequence of 
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where Zq is the canonical partition function, P = 1 /T is 
the inverse of the temperature, and k > is Grri^ times 
a geometrical number independent of Rq if Rq is large. 
The above inequality shows that the free energy grows 
with TV at least as N^^"^, and therefore is not extensive. 

The authors of Q argue that since to derive inequal- 
ity ^ I have introduced another length Rq, the relevant 
variable is 77 = Gm?N /TRq, and, since i?o ~ N^^^, we 
have T] ^ r]c, so that the system is deep in the collapse 
phase, where the results of do not apply. 

Obviously, this is wrong. Rq is not a characteristic 
length of the system. It is an auxiliary mathematical 
length, introduced just to prove that the diluted limit is 
ill-behaved. It has no physical meaning, and hence it is 
left unspecified. It can have any value, the only retriction 
being that it must scale with the number of particles 
as Rq ^ N^^^. By definition, the length entering 77 in 
equation (1) of ref. Q is the linear size of the spatial 
region available to the system. Hence, 77 = Gm? /TR. 
Only in such case the integrals over the coordinates that 
appear m equations (2.9), (2.15), etc., of |1| can be taken 
between and 1. Inequality Q can be written as 
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where k > is now a dimensionless purely geometri- 
cal number independent of the size of the system. Since 
R/Rq ~ iV^/'^, we see that the free energy scales N^^^ if 
77 is kept constant, and therefore is nonextensive. 

It is clear that the grand canonical partition function 
cannot scale as exp[Ng(ri, ^)] if the canonical partition 
function grows with N faster than exp[A^/(77)]. Hence, 
the grand canonical ensemble cannot give an extensive 
thermodynamic potential, either. 

Inequalities l(2Jl or (j2Jl are valid whatever the value of 
the temperature (i.e., of 77). Hence, they prove that the 
gas phase obtained in refs. 0, 01 does not exist in the 
diluted limit. For a finite system, we expect a gas phase 
at high temperature (small 77) and a collapse phase at low 
temperature (large rj), with the two phases separated by 
a phase transition or crossover at some rjc- As the system 
size grows, ?7c will decrease towards zero, the gas phase 
will shrink and the collapse phase will eventually cover 
the whole phase diagram. 
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In |3] it is shown that the diluted hmit exists order by 
order in the cluster expansion (a similar proof was given 
in for the high temperature expansion). This contra- 
dicts the proof of nonexistence of the diluted limit. Since 
inequalities (0) or ^ are derived rigorously, with no as- 
sumption, the high temperature and cluster expansions 
cannot be valid. The reasons why this kind of expansions 
fail have been analyzed in 2]. Basically, there are two 
possibilities: 

i) The series in r]/N may not converge in the N —^ oo 
limit, due to the contribution of high order diagrams that 
are naively suppressed by powers of 1 /N, but which actu- 
ally may give a significant contribution ought to the short 
distance divergences (the cut-off behaves as a = A/N, 
where A is a fixed length, Cf. Eq. (2.8) of the paper). 
Concerning this point, the statement that appears at the 
end of the introduction of the paper [sf, "one can take 
the limit N ^ oo and then a 0" does probably not 
hold due to the singularities of the integrals that give 
the coefficients of the cluster expansion. The rigorous 
limit is, obviously N —>■ oo, a — > 0, with Na fixed. The 
modification of the procedure to take the iV — ^ cxd limit 
may be at the core of the failure of the cluster expansion 
developed in 

ii) Even if the cluster expansion were convergent, the 
series could not represent the thermodynamic potential, 
since it is in contradiction with the rigorous result of 
0. In deriving the cluster expansion there is at least 
one mathematically unjustified exchange of limits that 
may invalidate the equality between the thermodynamic 
potential and the cluster series. The cluster expansion (I 



use the notation of Q), 

Qn{v) = ^ + Y.jf'^^+llj Z'^/'^' + ■ ■ ■ ' (4) 

is rigorous, since the number of terms in the sums is fi- 
nite for finite iV. To proceed further, the authors of Q 
expand /y in powers of rj/N , and exchange the sum and 
integral. Mathematically, it can be very difficult to ana- 
lyze the conditions under which this exchange of limits is 
allowed, but we can get some insight from physical intu- 
ition. Inequality © suggests that the system is collapsed 
for large iV if 77 is of order I respect to 1/N . This means 
that the cluster expansion is dominated by the higher 
order terms (many particles within a cluster). Hence, 
the canonical (Gibbs) integration measure on configura- 
tion space is very concentrated on collapsed configura- 
tions. Hence, this measure is very different from the flat 
measure d^ri that correspond to a gas. Keeping only 
the dominant of the expansion of fij in powers of rj/N 
means that one is using effectively the flat measure. To 
recover something similar to the concentrated true mea- 
sure, one has to sum an infinite number of terms in tj/N. 
In other words, the expansion in ry/TV is valid only in the 
gas phase. Inequality © suggests that the gas phase can 
only take place for of the order of N~'^/'^. Hence, the 
radius of convergence of the expansion in 77 shrinks to 
zero as iV — + 00. 

Similar statements claiming the existence of the di- 
luted thermodynamic limit have been made in 01 without 
even mentioning the results of ref. 2]. 
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